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ABSTRACT 

Heavy particles in galaxy clusters tend to be more centrally concentrated than light ones according 
to the Boltzmann distribution. An estimate of the drift velocity suggests that it is possible that the 
helium nuclei may have entirely or partially sedimented into the cluster core within the Hubble time. 
We demonstrate the scenario using the Navarro-Frenk- White profile as the dark matter distribution of 
clusters and assuming that the intracluster gas is isothermal and in hydrostatic equilibrium. We find that 
a greater fraction of baryonic matter is distributed at small radii than at large radii, which challenges the 
prevailing claim that the baryon fraction increases monotonically with cluster radius. It shows that the 
conventional mass estimate using X-ray measurements of intracluster gas along with a constant mean 
molecular weight may have underestimated the total cluster mass by ~ 20%, which in turn leads to an 
overestimate of the total baryon fraction by the same percentage. Additionally, it is pointed out that the 
sedimentation of helium nuclei toward cluster cores may at least partially account for the sharp peaks 
in the central X-ray emissions observed in some clusters. 

Subject headings: cosmology: theory — dark matter — galaxies: clusters: general — X-rays: galaxies 



1. INTRODUCTION 

An accurate determination of the total gravitating mass 
of galaxy clusters is crucial for the 'direct' measurements of 
the cosmic mass density parameter J7m through the mass- 
to-light technique (e.g. Bahcall, Lubin & Dorman 1995) 
and the baryon fraction method (e.g. White et al. 1993). 
The latter has received much attention in recent years be- 
cause of the rapid progress in X-ray astronomy and par- 
ticularly the large spatial extension of the hot and diffuse 
X-ray emitting gas in clusters. In the conventional treat- 
ment, the volume-averaged baryon fraction fb, defined as 
the ratio of the gas mass M gas to the total mass M tot of a 
cluster, is obtained by assuming a thermal bremsstrahlung 
emission and hydrostatic equilibrium for the intracluster 
gas. While such an exercise has been done for almost ev- 
ery X-ray selected cluster with known temperature, the 
resultant baryon fractions show rather a large dispersion 
among different clusters. In particular, fb appears to be 
a monotonically increasing function of cluster radius, and 
thereby cannot, in principle, be used for the determina- 
tion of the cosmic density parameter £Im since the asymp- 
totic form of fb at large radii does not approach a uni- 
versal value (White & Fabian 1995; David et al. 1995; 
Markevitch & Vikhlinin 1997a,b; White et al. 1997; Et- 
tori & Fabian 1999; Nevalainen et al. 1999; Markevitch 
et al. 1999; Wu & Xue 1999). This leads to an uncom- 
fortable situation that a greater fraction of dark matter is 
distributed at small scales than at large scales. It is com- 
monly believed that these puzzles have arisen from our 
poor understanding of the local dynamics inside clusters 
such as cooling/non-cooling flows, substructures, mergers, 
etc. Yet, a satisfactory explanation has not been achieved. 

The conventional cluster mass estimate from X-ray mea- 
surement of intracluster gas assumes a constant mean 
molecular weight /i. That is, all the chemical elements 



of the X-ray emitting gas obey exactly the same spatial 
distribution. Under this assumption, the total mass of 
the cluster M to t oc while the mass in gas M gas oc [i. 
Consequently the cluster baryon fraction /&(~ M gas / 'M to t) 
depends sensitively on the value of /i. It appears that our 
estimates of total mass and baryon fraction of a cluster 
would be seriously affected if its chemical elements do not 
share the same spatial distribution. 

The Boltzmann distribution of particles in a gravita- 
tional field cf> follows n oc exp {—m<j)/kT) 1 where n and m 
are the number density and the individual mass of the par- 
ticles, respectively. The heavier a particle is, the slower its 
thermal motion will be. So, heavy ions in the intraclus- 
ter gas will have a tendency to drift toward the cluster 
center. As a consequence of this sedimentation of heavy 
ions toward the central region, given that the cluster has 
survived for a sufficiently long time in the Universe, /i will 
no longer be a constant in the cluster. If the intraclus- 
ter gas is entirely composed of hydrogen and helium with 
their primordial abundances, the value of fi will be higher 
at cluster center while asymptotically reaches 0.5 at large 
radii for a fully ionized hydrogen gas. 

The chemically inhomogeneous distribution in galaxy 
clusters was initiated by Fabian & Pringle (1977), who 
studied the sedimentation of iron nuclei. They found that 
the iron nuclei in the X-ray gas may settle into the cluster 
core within a Hubble time. Taking into account the col- 
lisions of iron nuclei with helium nuclei, Rephaeli (1978) 
argued that a much longer sedimentation time than the 
Hubble time was required so that the iron nuclei would not 
have enough time to sediment into the core. In this Letter 
we wish to focus on helium nuclei instead of iron nuclei. 
Since the drift velocity vd oc AZ~ 2 (e.g. Fabian & Pringle 
1977), where A and Z are the atomic weight and the charge 
of the ion, respectively, the helium nuclei will settle much 
faster than the iron nuclei, and thus the sedimentation 
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may eventually take place. Indeed, Abramopoulos, Chana 
& Ku (1981) have calculated the equilibrium distribution 
of the elements in the Coma cluster, assuming an analytic 
King potential, and found that helium and other heavy 
elements are strongly concentrated to the cluster core. In 
particular, Gilfanov & Sunyaev (1984) have demonstrated 
that the diffusion of elements in the X-ray gas may signif- 
icantly enhance the deuterium, helium, and lithium abun- 
dances in the core regions of rich clusters of galaxies. For 
simplicity, we assume that the intracluster gas consists of 
hydrogen and helium, and then demonstrate their equi- 
librium distributions in the gravitational potential char- 
acterized by the universal density profile (Navarro, Frcnk 
& White 1995, 1997; hereafter NFW). As will be shown, 
under the scenario of helium sedimentation, the baryonic 
matter in clusters can be more centrally distributed than 
the dark matter, in contradiction to what is commonly be- 
lieved. This may open a possibility to solve the puzzle for 
an increasing baryon fraction with cluster radius. 

2. SEDIMENTATION OF HELIUM NUCLEI IN GALAXY 
CLUSTERS 

Our working model is based on the assumption that the 
intracluster gas composed of hydrogen and helium is in 
hydrostatic equilibrium at a temperature T with the un- 
derlying gravitational potential. This can be justified by a 
simple estimate of their relaxation time t eq (Spitzer 1978) 

im p m a k 3 ' 2 ( T T \ 3/2 

tcq ^ 8{2n) 1 / 2 n p Z 2 Z 2 e i \nA \r^ p ~ + ^J 

where m p , m a , Z p , and Z a are the masses and the charges 
of proton and helium nucleus, respectively. n p is the pro- 
ton number density, and the Coulomb logarithm is taken 
to be In A w 40. Obviously, t cq is much shorter than the 
present age of the Universe, and hence, the protons and 
the helium nuclei are readily in hydrostatic equilibrium at 
the same temperature. 

According to the Boltzmann distribution of particles in 
a gravitational field, heavy particles tend to be more cen- 
trally distributed in galaxy clusters than light particles. 
Therefore, with respect to protons, helium nuclei will tend 
to drift toward the cluster center. An immediate question 
is: Is the drift velocity sufficiently large for the helium nu- 
clei to have settled into the cluster core within the Hubble 
time ? The drift velocity vp of helium nuclei can be esti- 
mated through (Fabian & Pringle 1977) 

_ iA a ml ,2 g{2kTf/ 2 

V °~ 16^/2^2 e 4 nplnA 

= 8.8 x 10 6 cm/s f - a A 

1 \'i x 10- 8 cms^V 

X (l0-3cm-0 {wk) ' (2) 

where g is the gravitational acceleration. Eq.(2) indicates 
that within the Hubble time the helium nuclei can drift a 



distance 

r D = 1.8 Mpc ( 3 - 1Q j; cm g _ 2 ) ( IF ^ F s) \^ K \ 

(3) 

which is indeed comparable to cluster scales. Moreover, 
as t\d oc grip 1 , the value of t\d (and vd) increases rapidly 
with cluster radius r. Therefore, Eq.(3) suggests that the 
majority of the helium nuclei have probably sedimcntcd 
into the cluster core within the Hubble time. Note that, 
due to the requirement of electrical neutrality, the elec- 
trons of the same charge will simultaneously sediment 
along with the helium nuclei. Here, we have not consid- 
ered the effects of magnetic fields and subcluster mergers, 
which may somewhat retard the sedimentation of helium 
nuclei (Rephaeli 1978; Gilfanov & Sunyaev 1984). 

The sedimentation of helium nuclei in galaxy clusters 
will lead to a dramatic change of baryonic matter distri- 
bution. Consequently, the determination of gas and total 
mass of a cluster will be affected through the mean molec- 
ular weight ii. Another significant effect is that a sharp 
peak in the X-ray emission concentrated in cluster core will 
be expected due to the electron - helium nucleus radiation. 
This provides an alternative scenario for the 'cooling flows' 
seen in some clusters, for which a detailed investigation 
will be presented elsewhere. In the present Letter, we only 
focus on the dynamical effect. In the extreme case, the X- 
ray gas will be helium-dominated at cluster center, while 
hydrogen-dominated at large radii. The mean molecular 
weight n, which is commonly used as a constant, will be a 
decreasing function of the cluster radius. At cluster cen- 
ter, ii reaches 4/3, the value for a fully ionized helium gas, 
while at large radii, [i approaches 0.5, the value for a fully 
ionized hydrogen gas. 

The conventional mass estimate from X-ray measure- 
ment assumes a constant mean molecular weight of ti = 
0.59. As the total dynamical mass M to t of a cluster is 
uniquely determined by the intracluster gas at large radii, 
the difference between M t c ot (the value from the conven- 
tional method with \x — 0.59) and M tot as a result of 
helium sedimentation is simply 

M tot = ^M t c ot = 1.18M t c ot . (4) 

This indicates that the conventional method using X-ray 
measurement of intracluster gas, together with a constant 
mean molecular weight, may have underestimated the to- 
tal cluster mass by ~ 20%, which in turn, results in an 
overestimate of the total baryon fraction by the same per- 
centage. While the effect of the helium concentration to- 
ward cluster center can alter the gas distribution, the total 
mass in gas of the whole cluster remains unaffected because 
of the mass conservation. 

3. GAS DISTRIBUTION UNDER THE NFW POTENTIAL 

We now demonstrate how hydrogen and helium are dis- 
tributed in clusters described by the NFW profile 

9 {r/r s ){l + r/r s f {> 

As has been shown by Makino et al. (1998), such a poten- 
tial results in an analytic form of gas number density: 

n g as(x) - n gas (0) e-"-(l + z) %as/x , (6) 
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where x — r/r s , and 77 gas = 4irGnm p p s r^/kT. Except 
for the small core radius, n gas (x) is well approximated by 
the conventional (3 model. If we neglect the interaction 
between protons and helium nuclei, then 

n p = n p0 e- r >(l + x)^ x (7) 

n a = n QO e- 8 "/ 3 (l + x) 8 ^ (8) 

where r] = 2TrGm p p s r^/kT. Eqs (7) and (8) give n a oc 

rip , which differs from the prediction by Gilfanov & Sun- 
yaev (1984), n a oc n p . The discrepancy is probably due 
to the fact that we have not accounted for the diffusion- 
induced electric fields. 

We display in Fig.l the radial distributions of protons 
and helium nuclei as well as their combined result for a typ- 
ical nearby cluster with kT = 7 keV, ?7g as = 10 and r s = 1 
Mpc (e.g. Ettori & Fabian 1999). It is apparent that the 
intracluster gas is dominated by different elements at dif- 
ferent radius ranges: Within the core radius of a few tenth 
of r s , the number density of helium is about four times 
larger than that of hydrogen because of the sedimentation 
of helium nuclei, giving rise to a significant excess of both 
mass in gas (see Fig. 2) and X-ray emission in the central 
region of the cluster relative to the conventional model. 
Outside the core radius, helium profile shows a sharp drop 
and protons become the major component of the gas. It 
has been claimed that the total gas distribution can be ap- 
proximated by the (3 model (Makino et al. 1998). In fact, 
it is easy to show that both the narrow (helium) and ex- 
tended (hydrogen) components of intracluster gas can be 
fitted by the [3 models with different (3 parameters. This 
provides a natural explanation for the double (3 model ad- 
vocated recently for the cooling flow clusters (e.g. Ikebe 
et al. 1996; Xu et al. 1998; Mohr, Mathiescn & Evrard 
1999; etc.). Also plotted in Fig.l is the mean molecular 
weight calculated by p, — (n p + An a )/(2n p + 3n a ). The 
constant mean molecular weight is now replaced by a de- 
creasing function of radius. The asymptotic values of p at 
small and large radii are 4/3 and 1/2, respectively. 

We present in Fig. 2 a comparison of gas masses de- 
termined by the invariant and variant mean molecular 
weights using the same cluster parameters as in Fig.l. To 
facilitate the comparison, we require that the total par- 
ticles within r2oo remains unchanged, where r2oo is the 
radius within which the mean cluster mass density is 200 
times the critical mass density of the Universe (f2 = 1). It 
appears that the sedimentation of helium nuclei leads to a 



remarkable concentration of baryonic matter towards clus- 
ter center, which challenges the conventional prediction 
that the baryon fraction increases monotonically with ra- 
dius (Fig. 3). This opens a possibility that the asymptotic 
baryon fraction at large radii may match the universal 
value defined by the Big Bang Nucleosysthesis, although 
a detailed investigation will still be needed. 

4. CONCLUSIONS 

Intracluster gas is mainly composed of hydrogen and 
helium. Their average abundances over a whole cluster 
should be of the cosmic mixture. However, their spatial 
distributions are entirely determined by the underlying 
gravitational potential of the cluster, and thus follow the 
Boltzmann distribution. On the other hand, clusters are 
believed to have formed at redshift 2 ;> 1. Therefore, the 
helium nuclei in clusters may have entirely or partially 
scdimented into the central cores of clusters today. This 
will lead to a significant change of the radial distributions 
of gas and baryon fraction. In the present Letter, we have 
only discussed the impact on the dynamical aspect of clus- 
ters. We will present elsewhere the effect on the cluster 
X-ray emission (e.g., cooling flows, the double (3 model, 
etc.) as a result of the sedimentation of helium nuclei. 

In the conventional treatment where the mean molecular 
weight is assumed to be constant, one may have underesti- 
mated the total dynamical mass of clusters by ~ 20%. Us- 
ing a more vigorous way in which the NFW profile is taken 
as the background gravitational field of clusters, we have 
studied the hydrogen and helium distributions. Indeed, 
the sedimentation of helium nuclei toward cluster centers 
has significantly changed the distribution of intracluster 
gas, with gas being more centrally concentrated than dark 
matter. This may open a possibility to resolve the puz- 
zle that the baryon fraction increases monotonically with 
radius predicted by the conventional model. In a word, 
a number of cosmological applications of the dynamical 
properties of clusters will be affected by the sedimenta- 
tion of helium nuclei if it has really taken place during the 
evolution of clusters, which includes the determination of 
Qm through ft,, the constraints on the cosmological models 
through X-ray luminosity - temperature relation, etc. A 
detailed theoretical study, together with the observational 
constraints, will be made in subsequent work. 

We thank an anonymous referee for helpful comments. 
This work was supported by the National Science Foun- 
dation of China, under Grant No. 1972531. 



REFERENCES 



Abramopoulos, F., Chana, G. A., & Ku, W. H.-M. 1981 ApJ, 248, 
429 

Bahcall, N. A., Lubin, L., & Dorman, V. 1995, ApJ, 447, L81 
David, L. P., Jones, C, & Forman, W. 1995, ApJ, 445, 578 
Ettori, S., & Fabian, A. C. 1999, MNRAS, 305, 834 
Fabian, A. O, & Pringle, J. E. 1977, MNRAS, 181, 5P 
Gilfanov, M. R., & Sunyacv, R. A. 1984, Sov. Astron. Lett., 10, 137 
Ikebe, Y. et al. 1996, Nature, 379, 427 
Makino, N., Sasaki, S., & Suto, Y. 1998, ApJ, 497, 555 
Markevitch M., & Vikhlinin, A. 1997a, ApJ, 474, 84 
Markevitch M., & Vikhlinin, A. 1997b, ApJ, 491, 467 
Markevitch M., Vikhlinin, A., Forman, W. R., & Sarazin, C. L. 1999, 
ApJ, 527, 545 

Mohr, J. J., Mathiesen, B., & Evrard, A.E. 1999, ApJ, 517, 627 



Navarro, J. F., Frcnk, C. S., & White, S. D. M. 1995, MNRAS, 275, 
720 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 
Nevalainen, J., Markevitch, M., & Forman, W. 1999, ApJ, 526, 1 
Rephaeli, Y. 1978, ApJ, 225, 335 

Spitzer, L. Jr. 1978, Physics of Fully Ionized Gases (New York: 
IntsrscisncG J 

White, D. A., & Fabian, A. C. 1995, MNRAS, 273, 72 

White, D. A., Jones, C, & Forman, W. 1997, MNRAS, 292, 419 

White, S. D. M., Navarro, J. F., Evrard, A. E., & Frenk, C. S. 1993, 

Nature, 366, 429 
Wu, X.-P., & Xue, Y.-J. 1999, MNRAS, in press 
Xu, H., et al. 1998, ApJ, 500, 738 



4 



Fig. 1. — Distribution of intracluster gas tracing the gravitational potential by the NFW profile for a typical nearby cluster, kT = 7 keV, 
Vgas = 10 and r s = 1 Mpc. The saturated distributions of hydrogen and helium as well as their combined result are shown by dotted lines 
and solid line, respectively. All the profiles are scaled by the central number density of hydrogen. The corresponding mean molecular weight 
is displayed in the bottom panel by solid line. Also plotted is the conventionally adopted value fi = 0.59 (dotted line). 

Fig. 2. — Different mass components plotted against cluster radius for the same cluster in Fig.l. The dashed line shows the mass in gas 
obtained in the conventional model, in which the central gas density is assumed to be 10 — 3 cm -3 . The conservation of total baryonic particles 
(or mass) within a radius of r200 = 3.1 Mpc is required for both with and without the sedimentation of helium nuclei. 

Fig. 3. — A comparison of baryon fractions between the conventional model (dotted line) and the scenario with the sedimentation of helium 
nuclei (solid line). The cluster parameters are the same as in Fig. 2. 



